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Abstract 

Experimental evidence for the existence of strictly higher-order phase transitions 
(of order three or above in the Ehrenfest sense) is tenuous at best. However, there 
is no known physical reason why such transitions should not exist in nature. Here, 
higher-order transitions characterized by both discontinuities and divergences are 
analysed through the medium of partition function zeros. Properties of the distri- 
butions of zeros are derived, certain scaling relations are recovered, and new ones 
are presented. 



1 Introduction 



In its original format, the Ehrenfest classification scheme identifies the order of a phase 
transition as that of the lowest derivative of the Helmholtz free energy which displays 
a discontinuity there [T]. Typical transitions which fit to this scheme are first-order 
solid-liquid-vapour transitions and second-order superconducting transitions. There are, 
however, many transitions characterised by divergent rather than discontinuous behaviour. 
Examples include ferromagnetic transitions in metals and the spontaneous symmetry 
breaking of the Higgs field in particle physics, which display power-law or logarithmic 
divergent behavior as the transition is approached. The classification scheme has, in 
practice, been extended to encompass these scenarios and the order of a transition is 
commonly given by the order of the lowest derivative in which any type of non-analytic 
behaviour is manifest. 

It has long been suspected that transitions of Ehrenfest order greater than two (with a 
discontinuity at the transition point) do not exist in nature. However there is no obvious 
physical reason why this should be the case. In fact, recent experimental observations of 
the magnetic properties of a cubic superconductor have been ascribed to its possessing a 
fourth-order discontinuous transition (see also [3] where the existence of well defined 
anomalies in the specific heat at the transition point was claimed). A theory for higher- 
order transitions was developed in IH El E] and found to be consistent with experimental 
work. 

Higher-order phase transitions (with either a discontinuity or a divergence in an ap- 
propriate free-energy derivative) certainly exist in a number of theoretical models. There 
are third-order temperature-driven transitions in various ferromagnetic and antiferromag- 
netic spin models jTHH], as well as spin models coupled to quantum gravity [HI [101 ■ Recent 
theoretical studies also indicate the presence of third-order transitions in various supercon- 
ductors (111, DNA under mechanical strain |T21, spin glasses [12], lattice and continuum 
gauge theories and matrix models linked to supersymmetry [T^l. A fourth-order 
transition in a model of a branched polymer was studied in ^^1 and the Berezinskii- 
Kosterlitz-Thouless transition is of infinite-order [T7j . 

In this paper, we analyse higher-order transitions through the medium of partition 
function zeros. To set the notation, let t represent a generic reduced even variable and 
h be the odd equivalent so that t = T/Tc — 1 and h = H/kBT in the notation of the 
Ising model (i.e., T is the temperature, which is critical at Tc, and H is the external 
magnetic field). The critical point is given by {t,h) = (0,0). This may be the end-point 
of a line of first-order transitions, as is the case in the Ising or Potts models. In the 
Potts-like case where the locus of transitions is curved, we may instead assume that t 
and h are suitable mixed variables, so that h is orthogonal to t, which parameterizes 
arc length along the transition line JH]- The free energy in the thermodynamic limit is 
denoted by f{t, h) and its n*^-order even and odd derivatives are ft"'\t, h) and fh^\t, h) 
so that the internal energy, specific heat, magnetization and susceptibility are given (up 
to some inert factors) as e{t,h) = fl;^\t,h), C{t,h) = fl'^\t,h), m{t,h) = fj^\t,h), and 
x{t,h) = fj^\t,h), respectively. In the following, to simplify the notation, we drop the 
explicit functional dependency on a variable if it vanishes. 
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One then commonly describes as an m -order phase transition a situation where the 
first (m — 1) derivatives of the free energy with respect to the even (thermal) variable 
are continuous, but where the m**^ thermal derivative is singular, with a discontinuity or 
a divergence at the transition point. The lowest (m' — 1) derivatives of the free energy 
with respect to the odd (field) variable may also be continuous in t, with a singularity 
occuring in the m'*^ derivative. Thus a continuous specific heat is realized if m > 2 
and the susceptibility is also continuous if m' > 2 as well. This situation, which is not 
normally possible in a ferromagnet^, is the one analysed in [HE], in which m = m' > 2. 
Such higher-order transitions may be possible in branched polymers and diamagnets such 
as superconductors. 

In the more common scenario where m' is not necessarily the same as m, the scaling 
behaviour at the transition may be described by critical exponents at h = 0: 

ft\t)-r\ ft'\t)r~.t-^, fi'\t)^t^, (1.1) 

while, for the magnetization in field at t = 0, we write 

fi'\h)^h'/r (1.2) 

In the familiar case of a second-order transition (m = m' = 2), the exponents A and G 
become, in standard notation, a and 7, associated with specific heat and susceptibility, 
respectively. 

In the theoretical work of IS], the following scaling relations were derived for the case 
of a diverging higher-order transition in which m' = m: 

{m-l)A + mf3 + G = m{m-l) , (1.3) 
G = (3 {{m - 1)6 - 1) . (1.4) 

In the second-order case ()1.3|) and ()1.4|) become equivalent to the standard Rushbrooke 
and Griffiths scaling laws, 

a + 2/3 + 7 = 2, 7 = /3(5-l), (1.5) 

as one would expect. 

Since the seminal work by Lee and Yang [19j as well as by Fisher [20] , the analysis of 
zeros of the partition function has become fundamental to the study of phase transitions. 
Fisher zeros in the complex temperature plane pinch the real axis at the physical transition 
point. The locus of Lee- Yang zeros, in the complex magnetic-field plane, is controlled by 
the (real) temperature parameter and in the high-temperature phase, where there is no 

^ One can readily see this by considering the Rushbrooke scahng law (|1.5() together with hyperscaling 
which give "f/v = d — 2(3/v {d being dimensionality and v the correlation-length critical exponent). Since 
for a system of finite linear extent L, the magnetisation obeys (|m|) oc L~^/^ ^ and since completely 
uncorrelated ferromagnetic spins would lead by the central limit theorem to {\m\) cx L~'^l'^ , we obtain 
the bound ^jv < d/2, since the actual decay in the correlated case is slower. From this, one obtains the 
restriction that j/f cannot be negative for a ferromagnet. 
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transition, it ends at the Yang-Lee edge. Denoting the distance of the edge from the real 
axis by tyl, one has the generic behaviour 

rydt) ~ , (1.6) 

at a second-order transition. The exponent A is related to the other exponents through 
A = 275/(5-1). 

In the remainder of this paper, a number of results concerning the locus and density 
of zeros are presented. Higher-order transitions controlled by a single parameter are 
analysed in Sec. |21 where the locii and densities of the corresponding Fisher zeros are 
determined. Various restrictions on the properties of such transitions are established and 
simple quantitative methods for analysing them are suggested. In Sec. El the focus is on 
the zeros of the Lee- Yang variety where even and odd control parameters come into play. 
Here, the scaling relations (jl.Hj) and ()1.4|1 are recovered and elucidated and a number of 
other ones are presented. Finally, conclusions are drawn in Sec. EJ 



2 Fisher Zeros 

In the bulk of physical models the locus of Fisher zeros is linear in a suitable parameter, 
u, which is a function of t and can be parameterized near the transition point, Uc, by 

1211122112311211 

u{r) = Uc + r exp{i(j){r)) . (2.1) 

This singular line in the upper half-plane has < 0(r) < vr, while that in the lower half 
is its complex conjugate. 

In the thermodynamic limit, the (reduced) free energy is 

f(t) =2Re \n(u-u(r))g(r)dr , (2.2) 
Jo 

where g{r) is the density of zeros and is a cutoff. We are interested in the moments 
given by 

/i")(t) = 2i-ir-\n - l)!Re T r^^^dr , (2.3) 

JO [u — u{r)) 

and consider the cases of discontinuous and divergent m^'^-order temperature-driven tran- 
sitions separately. 

The difference in free energies on either side of the transition can be expanded as 
f+{t) — f-(t) = I]^iC„(u — Uc)^, where + and — refer to above and below Uc- For a 
discontinuous transition, c„ = for n < m, while 7^ and the discontinuity in the m^^ 
derivative of the free energy is 

A/f) = mlcm ■ (2.4) 

Now, the real parts of the free energies must match across the singular line (otherwise 
the transition would be of order zero) which, from ()2.1|) . means SJ^m ^n'""' cos?T,</)(r) = 0. 
Therefore the impact angle (in the upper half-plane), (p = limr^o 4>{r) , is 

+ / = 0,...,m-l. (2.5) 

2m 
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It is now clear that, under these conditions, vertical impact is allowed only at any dis- 
continuous transitions of odd order. A discontinuous second-order transition with impact 
angle 7r/2 is forbidden. Similarly an impact angle of 7r/6, for example, is only allowed at 
a transition of order 3 or 9 or 15, etc. This recovers disparate results for first-, second- 
and third-order transitions in ^U]; [2^] and ^01 which are associated with impact angles 
7r/2 (corresponding to / = 0), 7r/4 (/ = 0) and 7r/2 (/ = 1), respectively. The question 
now arises as to the mechanism by which the system selects its /-value. One expects that 
further studies of higher-order transitions will be required to provide an answer. 

Let t = u — Uc, T = te~^'^ and assume that the leading behaviour of the density of 
zeros is g{r) = gor^, where Qq is constant. If p is an integer, analytical extension of the 
integration (|2.3|) to the complex plane yields the following result for the n^^ derivative: 



P+i ^ 
ft\t)= -2in-l)\goReY.e-'''% 



(2.6) 



where 5 is a lower integral cutoff and 

_ p!rP+^-"ln(r - r) 

- {^n-l)\{p + l-n)\ ■ ^^-^^ 

One finds that all Tj terms vanish as the transition is approached from above or below, 
except the term for which j = p + 1. li n < p this term is constant and there is no 
discontinuity in f^^"^ across the transition, while for n = p + 1 it leads to a discontin- 
uous transition with A/^^^^^^ = 2TTgopl sin {p + 1)0. Therefore the first p derivatives are 
continuous across the transition while the {p + l)*'^ derivative is not. In other words, to 
generate a discontinous transition of order m under these assumptions, it is necessary and 
sufficient that p = m — 1, i.e., the leading behaviour of the density is 

9{r) = gor^-' . (2.9) 

From (j2.4|) and (j2.5j) . one now has = {—^Y'^'^go/fTi-, and the discontinuity in the rn"^ 
derivative of the free energy is related to the density of zeros as 

A/f) = (-l)'(m-l)!27r^?o. (2.10) 

This recovers the well known result that the latent heat or magnetization is related to the 
density of zeros at a first-order transition through Afj;^^ = 2ngQ |19j . 
We next consider an m^'^-order diverging transition where 

/f)(t)~|t|-^, (2.11) 

for 0<y4<l. Ify4 = 0, we are back to the discontinuous case or the case of a logarithmic 
as opposed to power-law divergence (see the discussion below), while if A > 1, it is more 
appropriate to consider the transition as (m — l)*'^ order. 
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Considerations similar to those in [231 may be used to show that in order to obtain 
the appropriate divergence it is necessary and sufficient that 



9{r) = gor'^-'-'' . (2.12) 

Indeed, from the general expression ()2.3|) . the form ()2.1ip is obtained provided (with 
r = tr') 

fR t^g(r)dr fR/n^-"'+^g(tr')dr' 
Re / , \ = Re / — — ./^ / — 2.13 
Jo {re't' - t)"^ Jo (r'e*<^ - 1)™ 

is independent of t as t ^ 0. The further condition that ^'(0) = gives A < m — 1. If 
m = 1, this violates the condition that < A < 1, leading to the requirement that m > 2 
for a diverging transition. On this basis, there are no diverging ffist-order transitions. 
This is consistent with experience. 

To demonstrate sufficiency, we put ()2.12|) into ()2.3p and use the substitution w = 
rexp {i(j))/\t\, to find, for the n**^ derivative of the free energy, 

ft\t) = goin - l)!|t|"-"-^e-^('"-^)'^/± , (2.14) 

in which 

rRe"l'/\t\ yjm-l-A 

/± = 2Re / -dw for t < . (2.15) 

Jo (l±w)" > 



If n < m, this vanishes as t ^ 0, establishing the continuity of the n^^ derivative there, 
while if n = m, one finds 

f^-\t) = -2g,\t\-^T{m-AnA)xh'''[^-'^lt .^ V.^n (2-16) 

In the case of a second-order transition, ()2.16|) recovers a result derived in ^Hl 1^ . Note 
that (j2.16p provides a direct relationship between the impact angle and the critical am- 
plitudes on either side of the transition. These critical amplitudes coincide if the impact 
angle is = {2N — A)tt /2{'m — A) where is any integer. In particular, if m is even an 
impact angle of 7r/2 results in the symmetry of amplitudes around the transition. This 
result was already observed in the second-order case in |23j. The implications of ()2.16|) are 
that, while this symmetry may be extended to all even-order diverging phase transitions, 
it does not hold for odd ones. 

If y4 = in ()2.12|) . the singular part of the m**^ derivative of the free energy becomes 

ft\t)=2im-l)\goxh;'^';'^l^\l. . , ^! J < ^ (2.17) 

^ 7 I (cos (m0) In |t| + TTsm (m0)) if t > . ^ ' 

This recovers a result in |231 if m = 2. Moreover, the discontinuity in the m*^ moment 
across the transition is consistent with ()2.5p and ()2.1Up . 

From ()2.9p and ()2.12|) . the integrated density of Fisher zeros is G{r) ~ f'^-^ (where 
^ = in the case of a discontinuous transition). For a finite system of linear extent L, 
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the integrated density is defined as Giitj) = (2j — 1)/2L'^ ^7\. Equating G(tj) to Giitj) 
leads to the scahng behaviour 

\tj\ ^ . (2.18) 

In the diverging case where hyperscahng (/(t) ~ C,{tY) holds, and m — A = 2 — a = ud, 
this recovers the usual expression, \tj\ L-^/", for finite -size scaling of Fisher zeros. In 
the discontinuous case, where A = 0, ()2.18|) yields 

This is a generalization of the usual formal identification of v with l/d, which applies to 
a first-order transition. Such a generalized identification was observed at the third-order 
(m = 3) discontinuous transition present in the spherical model in three dimensions [7j 
as well as in the Ising model on planar random graphs if the Hausdorff dimension is used 
for d [ini. 



3 Lee- Yang Zeros 

In the Lee- Yang case, where there is an edge, rYhif), to the distribution of zeros, the free 
energy is 

f{t,h) = 2Re ln{h- h{r,t))g{r,t)dr , (3.1) 

where the density of zeros is written as g{r,t) to display its t-dependency and where the 
locus of zeros is h{r,t) = rexp {i(j){r,t)). (If the Lee- Yang circle theorem holds, = tt/2 
and R = IT fTU].) The m'^^ field derivative of the free energy at = is 

fr\t) = 2i-ir'-\m' - iy.^^2lM$L 3±!^dx , (3.2) 

having used the substitution r = xryi^^t). As in the second-order case, we assume that 
^yl(''^) is sufficiently small near the transition point (t = 0) so that the upper integral 
limit diverges and compare with the limiting scaling behaviour in (jl.ip to find |221 

^7(r,t) = rVL(t)'"'-^$(^) , (3.3) 

where $ is an unknown function of its argument. Similar considerations yield, for the 
magnetization, 



ft\t, h) = 2t-«rYL(t)™'-iRe T -^^^dx , (3.4) 

rYL(t) 

which we may write as 

fl'\t,h) = t-^ry^tr'-'^,!^-^] . (3.5) 
6 



Comparison with ()1.2|) now gives "^{h/rYhit)) ~ (/i/'rYL(t)) • The t-dependence must 
cancel in ()3.5|) as t — » 0, giving the small-t scahng behaviour of the Yang-Lee edge under 
these circumstances to be 

as 

rYL(t) ~ t . (3.6) 

When m' = 2 and G = 7, this recovers the second-order transition behaviour of 
Furthermore, (j3.3j) now reads 

g{r, t) = tT^^^-<^(—^] , (3.7) 

and the expression for the magnetization in ()3.5p gives 

/«(t,/.) =t(^^^vl/^(^_^^ . (3.8) 

Strictly, this equation of state has been derived for t > 0, where there is an edge. However 
we may assume it can be analytically continued into the low temperature regime, where, 
taking the /i — ^ limit and comparing with the magnetization in , it yields the scaling 
relation 

In the situation where m' = this recovers the Griffiths- type scaling relation ()1.4p . 
derived in [Sj. 

Integrating 1)3.11) by parts gives, for the singular part of the free energy, 

where G{r, t) is the integrated density of zeros. From ()3.3|1 and ()3.6|) . the latter is G{r, t) = 
t^i^+^)/i'^'^'-^)^-^)F {r/rYhit)) in which F{x) = J^^{x')dx'. Again using r = xrYh{t) in 
()3.1())1 . and taking the upper integral limit to infinity, one has, for the free energy, 

f(t,h)=f^^r,[^^, (3.11) 

where T^{w) = 2Re F{x) / {we~^'^ — x)dx. The rn}^ temperature derivative of the zero- 
field free energy is therefore of the form //""V) ~ t^i^+^)/ii"^' -^)^-^)-^ . Comparison with 
(jl.lj) then yields the scaling relation 

A = m-G / + ^ . (3.12) 

[m — 1)0 — 1 

Together, ()3.9|) and ()3.12|) recover all four scaling relations derived in |3] in the more 
restrictive case where m' = m. In the second-order case (m = 2), they recover the 
standard Rushbrooke and Griffiths scaling laws of (|1.5p . 
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In fact, these laws also hold in the present case, albeit with negative a (and possibly 
7). To see this, let ft^^ it) ~ t""" and f^"* it) ~ t"^-^ (so that 0:2 = a and 72 = 7). Since 
ft"^\t) ~ t~^, one has, directly, that n — an = fri — A. Differentiating ()3.1H) with respect 
to field, now gives 

having used (j3.9|) and (j3.12j) and set h = 0. Now, one has 

7„ = (n - - /5, {n-l)an + n(3 + -fn = n{n-l), (3.14) 

which recover (jl.5|) when n = 2? 

The formute (jl.ip describe the behaviour of various moments as the critical point 
is approached tangential to the transition line (i.e., along h = Q). One may also be 
interested in the orthogonal behaviour, namely, the /i-dependence at t = 0. In the case of 
the /i-derivatives of free energy, this comes directly from ()1.2|) . For the t-derivatives, we 
may assume the power-law behaviour (at t = 0), 

fl'\h)^h^^ . (3.15) 

In the second-order case, ()3.15|) gives the /i-dependency of the internal energy and the 
specific heat at t = as e{h) = ft^\h) ~ h"" and C{h) = ft^\h) ~ These exponents 
are related to 5 and 7 through (see [21] and references therein) 

((^-1)(7 + 1) (^-l)(7 + 2) _ 
^ = 2 = J-^ 2. (3.16) 



Following the reasoning of |2T], we may argue that because there should be no phase tran- 
sition away from h = for any t, the free energy, f{t, h) in p.lip must be a power series in 
t there. So if J-'tf,{w) involves a term, w'^, the free energy involves -'<i)GS/((m -i)5~i)j^q 

which must be an integral power, A^, of t. This gives q = m' — [(m' — 1)5 — 1] (G + A^) / G5, 
or the power series 

/(t, h)=f: . (3.17) 

Differentiating appropriately, putting t = and comparing with ()3.15|) yields the scaling 
laws 

{m' -1)5- 1^ 
~G5 

In the second-order case with m' = 2 this recovers ()3.16|) with si = e and S2 = — o". 



=m'- ^'f \ G + j) . (3.18) 



^It is interesting to note the restrictions imposed on (5 at a higher-order transition coming from the 
first equation of H3.14|l . For n < m', 7„ should be negative, so, if (3 is positive, the best bound on S 
is S < \/{m! — 2). Also, the second formula in H3.14|l gives, for 2 < n < m — 1 and hence a„ < 0, 
6 > (rn — l)//3 — 1 or /3 > (m — l){m' — 2)/(m' — 1). These are no restraints in the familiar second-order 
case (where m = m' = 2 and large S is common), but are severe constraints at higher order. 
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4 Conclusions 



Different types of higher-order phase transitions have been analysed using the zeros of 
the partition function. In the Fisher case, the impact angle is restricted by the order and 
nature of the transition. For a transition with a discontinuity in ft™'\t), it is unclear 
how the system selects from the m permissible angles. For a divergent transition, the 
impact angle determines the relevant amplitude ratios. Finite-size scaling is seen to hold 
at higher-order transitions and the familiar formal identification of v with 1/d that is used 
at first-order transitions extends to u = m/d for discontinuous transitions of m**^ order. 

Lee- Yang zeros, on the other hand, are appropriate to the case where two parameters 
control the system. Here, they have been used as a route to derive scaling relations 
between associated even and odd exponents, which recover well known formulae in the 
second-order case, including the Rushbrooke and Griffiths laws. 

One of the main points of j2j is that many higher-order transitions may exist which 
have not yet been identified as such. Indeed determination of critical exponents or latent- 
heat-like discontinuities is notoriously difficult from numerical work on finite systems 
where there is no true transition and signals are smoothed out. There, amplitude ratios 
are often more discerning and here we see impact angles even more so, at least in theory. 
From the results herein, it would appear that analysis of the impact angle provides a very 
robust way to recognise the order of transitions. 



Acknowledgements: This work was partially supported by the EU RTN-Network 'EN- 
RAGE': Random Geometry and Random Matrices: From Quantum Gravity to Econo- 
physics under grant No. MRTN-CT-2004-005616. RK thanks Pradeep Kumar for an 
e-mail correspondence. 



References 

p. Ehrenfest, Commum. Kamerlingh Onnes Lab., Leiden, 20 Suppl. No. 75b (1933) 8. 

P. Kumar, Donavan Hall and R.G. Goodrich, Phys. Rev. Lett. 82 (1999) 4532. 

B.F. Woodfield, D.A. Wright, R.A. Fisher, N.E. Philips and H.Y. Tang, Phys. Rev. 
Lett. 83 (1999) 4622. 

P. Kumar, J. Low Temp. Phys. 106 (1997) 705. 

P. Kumar and A. Saxena, Phil. Mag. B 82 (2002) 1201; P. Kumar, Phys. Rev. B 68 
(2003) 064505. 

A.K. Farid, Y. Yu, A. Saxena and P. Kumar, Phys. Rev. B 71 (2005) 104509. 
T. Berlin and M. Kac, Phys. Rev. 86 (1952) 821. 



9 



[8] K. Kanaya and S. Kaya, Phys. Rev. D 51 (1995) 2404; H.G. Ballesteros, 
L.A. Fernandez, V. Martm-Mayor and A. Munoz Sudupe, Phys. Lett. B 378 (1996) 
207; Nucl. Phys. B 483 (1997) 707; A.P. Gottlob and M. Hasenbusch, J. Stat. Phys. 
77 (1994) 919; M. Kolesik and M. Suzuki, Physica A 216 (1995) 469; A.N. Rubtsov 
and T. Janssen, cond-mat / 0304447,. 

[9] V.A. Kazakov, Phys. Lett. A 119 (1986) 140; J. Ambj0rn, G. Thorleifsson and 
M. Wexler, NucL Phys. B 439 (1995) 187. 

[10] J. Ambj0rn, K.N. Anagnostopoulos and U. Magnea, Mod. Phys. Lett. A 12 
(1997) 1605; W. Janke, D.A. Johnston and M. Stathakopoulos, Nucl. Phys. B 614 
(2001) 494. 

[11] C. Cronstrom and M. Noga, Czech J. Phys. 51 (2001) 175; R. Werner, Phys. Rev. B 
67 (2003) 014505. 

[12] J. Rudnick and R. Bruinsma, Phys. Rev. E 65 (2002) 030902. 

[13] A. Crisanti, T. Rizzo and T. Temesvari, Eur. Phys. J. B 33 (2003) 203. 

[14] D.J. Gross and E. Witten, Phys. Rev. D 21 (1980) 446; M.R. Douglas and V.A. Kaza- 
kov, Phys. Lett. B 319 (1993) 219; M. Caselle, A. D'Adda, L. Magnea and S. Panzeri, 
Proc. 1993 Trieste Summer School in High Energy Physics and Cosmology (Trieste 
1993), ed. E. Cava et al, World Scientific, 1994; C.R. Gattringer, L.D. Paniak and 
G.W. Semenoff, Annals Phys. 256 (1997) 74; J. HaUin and D. Persson, Phys. Lett. 
B 429 (1998) 232; L. Li and Y. Meurice, |hep-lat/0507034[ 

[15] H. Fuji and S. Mizoguchi, Phys. Lett. B 578 (2004) 432. 

[16] P. Bialas and Z. Burda, Phys. Lett. B 384 (1996) 75. 

[17] V. Berezinskii, Zh. Eksp. Teor. Fiz. 59 (1970) 907 [Sov. Phys., JETP 32 (1971) 493]; 
Zh. Eksp. Teor. Fiz. 61 (1971) 1144 [Sov. Phys., JETP 34 (1971) 610]; J. Kosterlitz 
and D. Thouless, J. Phys. C 6 (1973) 1181; J. Kosterlitz, J. Phys. C 7 (1973) 1046. 

[18] F. Karsch and S. Stickan, Phys. Lett. B 488 (2000) 319. 

[19] C.N. Yang and T.D. Lee, Phys. Rev. 87 (1952) 404; ibid. 410. 

[20] M.E. Fisher, in: Lecture in Theoretical Physics, Vol. VIIC, ed. W.E. Brittin (Gordon 
and Breach, New York, 1968), p. 1. 

[21] R. Abe, Prog. Theor. Phys. 38 (1967) 72. 

[22] M. Suzuki, Prog. Theor. Phys. 38 (1967) 289; ibid. 1225. 

[23] R. Abe, Prog. Theor. Phys. 37 (1967) 1070; Prog. Theor. Phys. 38 (1967) 322. 
[24] M. Suzuki, Prog. Theor. Phys. 38 (1967) 1243. 



10 



[25] R.A. Blythe and M.R. Evans, Phys. Rev. Lett. 89 (2002) 080601; Braz. J. Phys. 33 
(2003) 464. 

[26] L.P. Widom, J. Chem. Phys. 43 (1965) 3892; ibid. 3898; L.P. Kadanoff, Physics 2 
(1966) 263; C. Domb and D.L. Hunter, Proc. Phys. Soc. 86 (1965) 1147. 

[27] W. Janke and R. Kenna, J. Stat. Phys. 102 (2001) 1211; W. Janke, D.A. Johnston 
and R. Kenna, Nucl. Phys. B 682 (2004) 618. 



11 



